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ABSTRACT 

By applying recent results for the slab correlation time scale onto cosmic ray scattering 
theory, we compute cosmic ray parallel mean free paths within the quasilinear limit. 
By employing these results onto charged particle transport in the solar system, we 
demonstrate that much larger parallel mean free paths can be obtained in comparison 
to previous results. A comparison with solar wind observations is also presented to 
show that the new theoretical results are much closer to the observations than the 
previous results. 

Key words: cosmic rays - turbulence - diffusion 



1 INTRODUCTION 

Cosmic rays (CRs) interacting with turbulent magnetic 
fields get scattered and accelerated (see Melrose 1968, 
Schlickeiser 2002). The theoretical description of these scat- 
tering and acceleration processes are essential for under- 
standing the penetration and modulation of low-energy cos- 
mic rays in the heliosphere, the confinement and escape of 
galactic cosmic rays from the Galaxy, and the efficiency of 
diffusive shock acceleration mechanisms. 

A key factor in CR scattering are the properties of the 
magnetic fields. A standard approach is the assumption of 
a superposition of a mean magnetic field -Bo ~ BqEz and 
a turbulent component 5B{x). Whereas the mean field can 
easily be meassured in the solar system (here we find approx- 
imately _Bo ~ 4 — 5nT), the turbulent component has to be 
emulated by turbulence models. In the literature there is 
no consensus available about the true turbulence properties 
(see Cho & Lazarian 2005 for a review). In the solar system, 
however, some turbulence properties such as the wave spec- 
trum can be obtained from meassurements (see e.g. Denskat 
& Neubauer 1983, Bruno & Carbone 2005). 

More unclear are the orientation of the turbulence wave 
vectors (also refered to as turbulence geometry) and the dy- 
namical decorrelation of the magnetic fields. In a recent CR 
diffusion study (Shalchi et al. 2006) a slab/2D composite 
model was combined with a nonlinear anisotropic dynami- 
cal turbulence (NADT) model. This model can be used to 
reproduce meassured CR mean free paths parallel and per- 
pendicular to the mean field Bq. The authors of this article 
assumed that the slab correlation time scale is independent 
of the wave vector k. 

In a recent study (Lazarian & Beresnyak 2006), how- 
ever, it was shown that the slab time scale is indeed 



fc— dependent More precisely it was found that ~ "(c = 
VA\fk\jL\ here we used the correlation time fc, the corre- 
lation rate 7^, the Alfven speed va, and the outer scale of 
the turbulence L. It is the purpose of this article to apply 
this new result of the slab correlation time scale onto cosmic 
ray parallel diffusion. A comparison with solar wind obser- 
vations of the parallel mean free path is also presented. It is 
demonstrated that we can find a much larger parallel mean 
free path if we employ the correlation time scale of Lazarian 
& Beresnyak (2006). 

In Section 2 we explain the turbulence model that is 
used in this article. In Section 3 a quasilinear description 
of cosmic ray scattering is combined with this turbulence 
model to derive analytic forms of the pitch-angle diffusion 
coefficient and the parallel mean free path. In Section 4 we 
evaluate these formulas numerically to compute diffusion co- 
efficients and we also provide a comparison with previous 
results and solar wind observations. In the closing Section 5 
our results are summerized. 



2 THE TURBULENCE MODEL 

The key input into a cosmic ray transport theory like QLT 
is the tensor P;,„ which describes the correlation of the tur- 
bulent magnetic fields; 



,(fc,t) = (5B,(fc,t)5B;;(fc,o)) 



(1) 



Therefore the fc-dependence and the time-dependence of the 
correlation tensor Pim{k,t) have to be specified which is 
done in the next two subsections. 

^ That study put to the test the idea of the damping of slab per- 
turbations by the ambient turbulence in Yan &i Lazarian (2002), 
Farmer & Goldreich (2004). 
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Energy range 



Inerlial range 



10 , 10 



Dissipation ranged 



10 10 k 10 



Table 1 . Previous models for the slab dynamical correlation func- 
tion. Here Vj^ is the Alfven velocity, a is a parameter which allows 
to adjust the strength of dynamical effects, and /3 = aZ/l2D is 
used in the NADT (nonlinear anisotropic dynamical turbulence)- 
model (Z = V25B2d/V^^ = V2vaSB2d / Bq). 



Model 



r(fe,i) 



Magnetostatic model 1 
Damping model of dynamical turbulence g,—°''"A\k\t 
Random sweeping model g— (au^fct) 



Plasma wave model for shear Alfven waves 
NADT-model 



-/3t 



Figure 1. The power spectrum used for our calculations. We used 
a general spectrum with energy, inertial, and dissipation range. 
The dissipation wavenumber divides the inertial range from 
the dissipation range. 



2.1 The slab model for the turbulence geometry 

For mathematical simplicity we employ the often used slab 
model for the turbulence geometry. Physical processes that 
can induce slab modes are differential damping of fast mode 
waves (Yan & Lazarian 2002) or instabilities (Lazarian & 
Beresnyak 2006). 

By assuming the same temporal behavior of all tensor 
components we have 

Pi^{k,t) = PiT{k)-T'"^\^,t) (2) 

with the dynamical correlation function r*'''''(fe|| , t) and the 
time independent correlation tensor Pi^^ik). The tensor 
Pim^{^) is determined by the turbulence geometry and the 
wave spectrum whereas the function r'''"'(fc, t) describes dy- 
namical effects. 

For pure slab turbulence we have 5Bi{x) = 5Bi{z) and 
therefore 

5(fc±) 



Slab / 7 

g {k\ 

Sim — 



ki kr 







if l,m=x,y 
if 1 or m = z 



(3) 



In Eq. ([3)l we assumed the case of vanishing magnetic helic- 
ity. The function 3'''"''(fc||) is the slab wave spectrum which 
can be approximated by a power-law spectrum with energy, 
inertial, and dissipation range (see Fig. 1): 

slab/, \ Civ) 2 

for s; < k 

(1 + k\lli^^y" for fc„,„ s: < fcd(4) 

(l + fcd'Li,)"''(fcdA||)'' for fcd^fcii 

where we used the normalization function 
1 V{v) 



C(v) 



2^^ r(i/ - 1/2) ■ 



(5) 



Furthermore, we used the slab bendover scale Isiah, the dis- 
sipation wavenumber kd, the turbulence strength SB^i^,,, the 
inertial range spectral index s = 2i/ and the dissipation 
range spectral index p. The parameter kmin = indi- 



cates the smallest wavenumber. Also the dynamical correla- 
tion function r°'"''(fe|| , t) has to be determined which is done 
in the next subsection. 



2.2 Improved form of the slab correlation time 
scale 

In the past several models have been developed to approx- 
imate the dynamical correlation function r'''"'(fc|| , t) (e.g. 
Schlickeiser & Achatz 1993, Bieber et al. 1994, Shalchi et al. 
2006). Some examples are given in Table [1] As in the most 
previous models we assume an exponential decorrelation of 
the magnetic fluctuations and, therefore 



In a recent article (Lazarian c 
that for slab modes 



and therefore 



r''''''(fc||,t) 



(6) 

: Beresnyak 2006) is was shown 

(7) 

(8) 



with the Alfven speed va, the plasma box size L, and with 
the plasma wave dispersion relation of shear Alfven waves 



iO = jVAk^\. 



(9) 



The form of Eq. ([7]) can be justified by the following ar- 
gument (for a more detailed explanation see Lazarian & 
Beresnyak 2006) ; consider a wavepacket of Alfven waves that 
moves nearly parallel to the magnetic field with the disper- 
sion of angles Skj_/k\^ ~ 0^. The individual waves follow the 
local direction of the magnetic field lines. As a result, the 
dispersion in angles of the wave packet cannot be less than 
the dispersion of angles due to the ambient Alfven turbu- 
lence, Qk > Obk- The latter for the Goldreich-Sridhar (1995) 
model 



,2/3.-1/3 



5v ~ VA{k±LA)^^^ 



(10) 



is Obfe ~ 5Bk/Bo ~ {k±L)~^^^ . The modes with minimal 
are the fastest growing ones. As we establish below (see Eq. 
(|11|) ') they are the least damped. Therefore for our simplified 
treatment we shall limit our attention to the wavepackets 
with resonant k'^^ ^ rp and Qbk ~ Qk- One can determine 
the characteristic perpendicular wavenumber k± ~ Sk± ~ 
{rp/ Ly^* of the most parallel modes that are created by 
streaming CR's. The strong Alfvenic turbulence decorrelates 
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the wavepacket with fcx on the time scale wxfcj_. Thus using 
the above expression for k± and Eq. (|10|) we get 

^ = 7c ~ -kl_V±_ ~ -VAk^l^ L""^''^ ~ -VATp'^^'^ L"'^''^ (11) 

and thus Eq. (O, which, up to the "-" that we used to denote 
the damping nature of the process, coincides with the damp- 
ing rate obtained in Farmer & Goldreich (2004) and with 
the results of the numerical simulations shown in Lazarian 
& Beresnyak (2006, Fig. 1 therein). 

The parameter j is used to track the wave direction 
(j = +1 is used for forward and j = — 1 for backward to 
the ambient magnetic field propagating Alfven waves). A 
lot of studies have addressed the direction of propagation of 
Alfvenic turbulence, see for instance Bavassano (2003). In 
general one would expect that closer to the Sun the most 
waves should propagate forward and far away from the Sun 
the wave intensities should be equal for both directions. In 
the current paper we are interested in turbulence parameters 
at 1 AU. Thus, we simply assume that all waves propagate 
forward and we therefore set j = +1. 

In the following we employ the dynamical correlation 
function defined by Eqs. ([8]) and ((9)1 to determine the parallel 
mean free path within the quasilinear approach. 



3 THE QUASILINEAR PARALLEL MEAN 
FREE PATH 

In the current paper we employ quasilinear theory (QLT, 
Jokipii 1966) to calculate the parallel mean free path. QLT 
can be seen as a first order perturbation theory in the small 
parameter SB/Bq. Whereas it was shown previously (e. g. 
Michalek & Ostrowski 1996) that QLT is accurate even if 
SB ~ Bo if we assume slab geometry and a wave spectrum 
without dissipation range, it was realized by more recent 
test particle simulations that for non-slab models and for 
steep wave spectra nonlinear effects are important and QLT 
is no longer accurate (see e.g. Qin et al. 2006 and Shalchi 
2007). In the current article we only investigate the slab 
model and we can therefore assume validity of QLT as it is 
shown in Shalchi et al. (2005) that the unphysical contri- 
bution of large scale motions arising in QLT due to its in- 
ability to account for the conservation of adiabatic invariant 
is small. According to Jokipii (1966), Hasselmann & Wib- 
berenz (1968), Earl (1974), and Shalchi (2006), the parallel 
mean free path results from the pitch-angle-cosine average 

as 



of the inverse pitch-angle Fokker-Planck coefficient D^^ 



All — — 



dfj. 



(12) 



The pitch-angle-cosine ^ is defined as = Accord- 
ing to Teufel & Schlickeiser (2002, Eq. 25) the pitch-angle 
Fokker-Planck coefficient can be written as 



»^(1-^" 
2BI 



i / k±v± 
1+1 I 

A' 

Jn+l 



V 


Q , 


k 






a } 


k_ 





PlUk) 



Jn-l 



k±v± 



X P^^(fc)e+^'-+PiB(fc)e 



(13) 



if we use helical coordinates 

SBl = {SB^ + iSBy) , 

v2 

SBr = -^{SB^^iSBy) 
v2 

and if we neglect electric fields 0. In Eq. (I13p we used the 
resonance function 



(14) 



-i{k\\v\\ +nn)t -pi 



r{k,t) 



(15) 



K(fc) ^ Re dt e 
Jo 

where the index i stands for the different turbulence mod- 
els. For pure slab geometry we have according to Teufel & 
Schlickeiser (2002, Eq. 33) 

PRR{k) = p^^(fc) =g="'''(fc||)^M, 



PRLik) = PLRik)=0 



(16) 



if we assume vanishing magnetic helicity. For the pitch-angle 
Fokker-Planck coefficient we then find 

r^3lab TvQ'^ (1 — fl^) f slabt, \ 

D„„ = ^^r:;, I g (k\^) 



^0 



(17) 



The resonance function for pure slab turbulence has the form 

poo 

Ri''\k«)=Rel dte-"''=ni+"">*-r''''(fc||,i). (18) 



With Eq. (O for r^'''''(fc|| , t), the integral in Eq. ^ is ele- 
mentary and we obtain 



Rt 



(19) 



With this Breit-Wigner-type resonance function the slab 
Fokker-Planck coefficient can be written as 



^0 



dk\\ g''^''^(k\\) 



VAyJk\\/L 



^ v\k«/L+{k«v«+nQ. — VAk\\ 

n — ±1 



(20) 



With the integral transformation x = lsiabk\\ and with the 
parameters R — Rh/lsiab ~ v/{^llsiab) and e — va/v we 
obtain 

slab _ 27r(l-/i^) 



D 



^0 tslab 



I 1 slab ( J ^ 

■ / dxg I fc|| = 

^ Q \ I' slab 



E 



VA\Jk\\IL 



„tii (VA^^JLI^Y + {xR{p. - e) + n\ 
The slab spectrum of Eq. Q can be written as 



(21) 



Because of the high conductivity of cosmic plasmas, there are 
no large-scale electric fields = Eq = and we thus have 

B = BoBz +SB, E = SE with the turbulent electric and magnetic 
fields {SE, SB). The reason for using the model of purely magnetic 
fluctuations is that the electric fields are much smaller than the 
magnetic fields, since we have SEi ~ —SBj <C SBj, with the 
Alfven speed va which is much smaller than the speed of light, c. 
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9 

with 



(22) 



7 Slab / J- \ 

h {X,Xmin,V,t,,p) 

for < a; < 

(1 + x^)-" for x^i„ ^x <^ (23) 

for 



where we used ^ = Isiabkd- Then we find for the dimension- 
less Fokker-Planck coefficient Dtf = D^'^*" 



MM 

,2^ xd2 



f^slab ^ C{v){l-ll')SBii,, 



da; h^^'^^{x,Xmin,J^,£,,v) 



E 



UA y^x/{lslabL)/ii 



„t^i (^Av/^/(/.ia6L)/f^)2 + [a;7?(/i - e) + n] = 
By using 

VA 



VA V 



Islab , 

5' V IslabJ-J V S'tslab V -tJ 



(24) 



(25) 



we finally find 

C(i.)(l-M^)5BL. 



^0 



da; h''^°'^{x,Xmin, v,^,p) 



eRy/xrj 



(ei?V5^)2 + [xi?(M-e)+n]^ 



(26) 



In Eqs. (|22|l - (|26|) we used p = hiab/ho and a;min = V = 
Isiab/L. A numerical investigation of the integral of Eq. (f26] 
is presented in the next section. 



4 NUMERICAL RESULTS FOR PITCH-ANGLE 
AND PARALLEL DIFFUSION 

In this section we evaluate the formulas for pitch-angle diffu- 
sion (Eq. (|26[) ) and for the mean free path (Eq. (|12p ') derived 
in the last sections numerically for the parameter-set of Ta- 
ble [5] which should be appropriate for interplanetary condi- 
tions at 1 AU heliocentric distance. All formulas depend on 
the parameter e which can be expressed as (see Shalchi et 
al. 2006) 



VA 
V 



VA 



R 



with 



Istab ■ Bo 



0.511MV for e" 
938MV for 



(27) 



(28) 



For the heliospheric parameters considered in the current 
paper we have for electrons Ro{e~) ~ 9.2 • 10~^ and protons 
-Ro(p+) ~ 0.169. 

In order to determine the transport coefficients of the 
isotropic part of the particle distribution function, e.g. the 
parallel mean free path, we must restrict our calculations to 



Table 2. Parameters used for our numerical calculations. The 
values should bo appropriate for the solar wind at 1 AU helio- 
centric distance. If a parameter is different from the values below 
we note this separately in the corresponding figures and discus- 
sions. These values correspond to the values used in Shalchi ct al. 
(2006). 



Parameter 


Symbol 


Value 


Inertial range spectral index 


2v 


5/3 


Dissipation range spectral index 


P 


3 


Alfven speed 


HA 


33.5 km/s 


Slab bendover scale 


^slab 


0.030 AU 


Slab dissipation wavcnumber 


^slab 


3- 10« (A(7)-i 


Mean field 


Bo 


4.12 nT 


Turbulence strength 


^Bslab/Bo 


1 




Figure 2. The pitch-angle Fokker-Planck coefficient for protons 
(dashed line) and electrons (solid line). Also the magnetostatic 
dissipationless model is shown (dotted line). Visualized is a lin- 
ear plot of the dimensionless pitch-angle Fokker-Planck coefficient 
Df^fj, as a function of the pitch-angle cosine fi. All results are for 
pure slab geometry and for small rigidities R = 10~*. 



t = va/v ^ 1 (see Schlickeiser 2002 for a detailed explana- 
tion). Thus, we can only consider rigidities which satisfy the 
following condition: 



Ro 



Ro 



VA 



(29) 



VA) 



For VA ~ 33.5 km/s we find for electrons the restriction 
R{e~) > 10"** and for protons i?(p+) > 2 • 10"^. 



4.1 The pitch-angle Fokker-Planck coefficient D^jp 

Fig. [2] shows the (dimensionless) pitch-angle Fokker-Planck 
coefficient 1),^^ = D^^fj^lsiab/v as a function of the pitch- 
angle-cosine fi. For the rigidity we assumed R = 10~ . In this 
case the parameter e (see Eq. Ip7)l ') is approximately e ^ 0.2. 
It seems that the minimum of the pitch-angle Fokker-Planck 
coefficient for protons can be found at /i ~ e. In general the 
pitch- angle Fokker-Planck coefficient is no longer equal to 
zero at 90° (/x = 0) as in the magnetostatic model, so that we 
no longer obtain an infinitely large parallel mean free path 
as in magnetostatic models. It should be noted, however. 
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lO"'^ 10"" 10"^ 10""^ 10' 10° 10^ 



L slab 

Figure 3. The quasilinear parallel mean free path Ay for x^in = 
77 = 0.1 for electrons (solid line) and protons (dashed line), and 
for Xmin = rj = 0.01 for electrons (dash-dotted line) and pro- 
tons (dotted line) in comparison with different solar wind obser- 
vations: Palmer consensus (Palmer 1982, box), Ulysses observa- 
tions (Gloeckler et al. 1995, dot), AMPTE spacecraft observations 
(Mobius et al. 1998, vertical line). 

that QLT itself is questionable close to 90°. By considering 
Fig. [2] we find that at least for protons pitch-angle scatter- 
ing close to 90° is very strong due to the dynamical effects. 
Therefore one could assume that nonlinear effects which also 
lead to nonvanishing pitch-angle scattering at 90° could be 
neglected. 

4.2 The parallel mean free paths A 

Here we present theoretical results for the parallel mean free 
path. We compare our results with the Palmer consensus 
(Palmer 1982) and pickup ion observations (Gloeckler et al. 
1995, Mobius et al. 1998). In Fig. |3] we show the parallel 
mean free path in comparison with observations. We com- 
puted the parallel mean free paths for two different values 
of the smallest wave number Xmin, namely for Xmin = 0.1 
and Xmin = 0.01. In addition to the Palmer (1982) results 
we compare our results also with pickup ion observations: 

(i) Gloeckler et al. (1995) concluded from Ulysses obser- 
vations that the parallel mean free paths of pickup protons 
is 2 AU at 2.4 MV rigidity (they stated conservatively that 
A II is of order 1 AU but actually they obtained the best fit 
for 2 AU). It should be noted that this observation was at 
high heliographic latitudes, and at a heliocentric distance 
of 2.34 AU; these differences should be remembered when 
comparing with observations at Earth orbit. 

(u) Mobius et al. (1998) concluded from AMPTE (Active 
Magnetospheric Particle Tracer Explorers) spacecraft obser- 
vations that the parallel mean free paths of pickup helium 
ranges from 0.16 to 0.76 AU at 5.6 MV rigidity in the data 
they analyzed. 

Both results are also illustrated in Fig. [3] As shown we can 
reproduce the observations theoretically. 




L slab 

Figure 4. The parallel mean free path Ay for Xmin = t] = 0.1. 
Shown are QLT results for electrons (solid line) and protons 
(dashed line) in comparison with the theoretical results obtained 
by Shalchi et al. (2006) by applying the NADT-model in combina- 
tion with the slab/2D model for the turbulence geometry (results 
for electron are represented by the dotted line, and the results for 
protons are represented by the dash-dotted line). 

In Fig. 2] we compare our new theoretical results with 
the results obtained by Shalchi et al. (2006) by applying 
the NADT-model in combination with the slab/2D model 
for the turbulence geometry. As shown, for small rigidities, 
where dissipation effects are important, we obtain a much 
larger parallel mean free path for electrons. For medium 
rigidities where the charged particles interact resonantly 
with the inertial range of the spectrum, the new parallel 
mean free path is about a factor 5 smallei[3 since we as- 
sumed pure slab fiuctuations (in comparison to the slab/2EQ 
model used in Shalchi et al. 2006). It seems that the im- 
proved model for the slab correlation function provides a 
much larger electron parallel mean free path for low ener- 
getic particles. 



5 SUMMARY AND CONCLUSION 

The theoretical explanation of measured parallel mean free 
paths in the solar system is a fundamental problem of space 
science. In a recent article (Shalchi et al. 2006) it has been 
demonstrated the these observations can indeed be repro- 
duced theoretically. By using recent results of turbulence 
theory (Lazarian & Beresnyak 2006) we further improved 

That fact that the parallel mean free path is much smaller for 
pure slab geometry (in comparison to the 20% slab / 80% 2D 
result of Shalchi et al. 2006) and medium rigidities, can be seen 
as a trivial result, since within QLT we find no gyro-resonant 
scattering due to 2D modes for magnetostatic turbulence. 
■* The slab/2D model used in Shalchi et al. 2006 assumes a super- 
position of pure slab and pure 2D modes. It is well known that this 
model can only be approximatelly correct, since real turbulence 
has a certain distribution of wave vectors in the fc— space. 
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the dynamical correlation function which is a key input in 
transport theory considerations. It is demonstrated in this 
article that the improved slab correlation time scale (see Eq. 
0) leads to a much larger parallel mean free path (see Fig. 
[3|. This effect is important since it was argued in several pre- 
vious articles that the theoretical parallel mean free path is 
too small (Palmer 1982, Bieber et al. 1994) in comparison 
with solar wind observations. 

Another problem of cosmic ray scattering theory is the 
importance of nonlinear effects. Whereas we have applied 
QLT in the current article it was argued in other papers 
(e.g. Shalchi et al. 2004) that nonlinear effects are impor- 
tant for parallel diffusion. However, these nonlinear effects 
are directly related to the interaction between charged parti- 
cles and 2D modes. These 2D modes were neglected since we 
assumed pure slab fluctuations. Therefore, QLT can be ap- 
plied and the results presented in this article should be valid. 
For non-slab models, where 2D modes are present, however, 
the applicability of QLT is questionable. It has to be sub- 
ject of future work to explore the validity of QLT for realistic 
turbulence models such as dynamical turbulence models in 
non-slab geometry. 
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